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Abstract. This article is concerned with the stabihty and long- 
time dynamics of structures arising from a structureless state. The 
paradigm is suggested by developmental biology, where morpho- 
genesis is thought to result from a competition between chemical 
reactions and spatial diffusion. A system of two reaction-diffusion 
equations for the concentrations of two morphogens is reduced to 
a finite system of ordinary differential equations. The stability of 
bifurcated solutions of this system is analyzed, and the long-time 
asymptotic behavior of the bifurcated solutions is established rig- 
orously. The Schnakenberg and Gierer-Meinhardt equations are 
discussed as examples. 



1. Introduction 

Morphogenesis — that is, the development and formation of tissues 
and organs — is one of the main mysteries in hving organisms. How 
does structure emerge from a structureless state without the apparent 
action of an external organizing force? A major factor seems to be 
the competition between chemical reactions and spatial diffusion of 
substances called morphogens, which are present in the cells. The idea 
goes back to the pioneering work of Turing in 1952 [23], who noted 
that diffusion in a mixture of chemically reacting morphogens can cause 
instability of a spatially uniform steady state and lead to the formation 
of spatial patterns; see the article by Cross and Hohenberg [3] and the 
recent text by Hoyle |8j for a comprehensive overview of the theory of 
pattern formation and Turing analysis. 

Turing's analysis, which is essentially a linear eigenvalue analysis, 
has been a basic tool in the study of nonlinear reaction-diffusion sys- 
tems; in fact, it has provided insight into the behavior of nonlinear 
systems as well, since the latter can often be approximated, at least 
for brief lengths of time, by linearized systems. But as time evolves, 
the nonlinear structure takes over, and other tools are needed to study 
the long-time behavior. In certain cases, where the existence of in- 
variant regions for reaction-diffusion systems can be established, the 
solution of a nonlinear system remains bounded [22]; however, even in 
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these cases it is not known rigorously whether the patterns persist in 
the long run, even though the idea is supported by many numerical 
simulations. 

The mathematical literature contains many instances of weakly non- 
linear stability analyses for reaction-diffusion systems. An early refer- 
ence is [6] , where a center-manifold approach is used; other, more recent 
references are [26l [27]. Sometimes, special techniques have been ap- 
plied to the study of Turing patterns in different regimes. For example, 
Ref. |n] deals with the stability of symmetric iV-peaked steady states for 
systems where the inhibitor diffuses much more rapidly than the activa- 
tor. We also mention Refs. [H [12] , which deal with the Schnakenberg 
model on a two-dimensional square domain, where spatially varying 
diffusion coefficients cause the removal of the degeneracy of the Turing 
bifurcation. 

Weakly nonlinear stability analyses can be justified rigorously on 
the basis of modulation theory and a Ginzburg-Landau approximation; 
see, for example, Refs. [21 [H [201 [211 [2S] ■ Murray's monograph [H] gives 
applications to biological systems such as animal coat patterns. 

The purpose of the present work is to investigate the nonlinear sta- 
bility and persistence of spatiotemporal patterns on bounded domains. 
The investigation is based on recent results of Ma and Wang [11] on 
attractor bifurcation for nonlinear equations. The attractor bifurcation 
theorem (see the Appendix, Section IA.1|) sums up the basic features 
of a stability-breaking bifurcation. Starting from the original partial 
differential equation, it identifies and characterizes the local basins of 
attraction based on the multiplicity of the eigenvalues near a bifurca- 
tion point. Thus, the attractor bifurcation theorem gives the complete 
picture, rather than the caricature given by the amplitude equations. 

A second essential feature of the present investigation is a center- 
manifold reduction to reduce the partial differential equation to a finite- 
dimensional dynamical system. The reduction requires the computa- 
tion of the center-manifold function and the interaction of the higher- 
order eigenfunctions with the eigenspace belonging to the leading eigen- 
values. Such a reduction is inherently difficult, and for this reason one 
usually resorts to a generic form of the reduced equation which is some- 
what detached from the original. On the other hand, a center-manifold 
reduction offers a practical way to find the structure of the local at- 
tractors of the original partial differential equation. These attractors 
completely describe the local transitions, and their basins of attraction 
define the long-time dynamics associated with the transitions. Since 
these are exactly the features of interest, we have taken this approach 
and focused much of our efforts on the center-manifold reduction. As 
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a result, we are able to characterize the types of transitions in terms 
of explicitly computable parameters which depend only on the domain 
and the values of the physical parameters of the system under consid- 
eration. 

We prove that spatiotemporal patterns in reaction-diffusion systems 
of the attractor-inhibitor type can arise as the result of a supercritical 
(pitchfork) or subcritical bifurcation. The former results in a contin- 
uous transition, the latter in a discontinuous transition. In the case 
of diffusion on a (bounded) interval or on a rectangular (non-square) 
domain, we prove that the attractor consists of two points, each with 
its basin of attraction (Theorem 14. Fig. [2]). In the case of diffusion 
on a (bounded) square, the phase diagram after bifurcation consists 
of eight steady-state solutions and their connecting heteroclinic orbits 
(Theorem 15. ![ Fig. [3]). The conditions for the stability of these bifur- 
cated steady states and the heteroclinic orbits are explicit; they can 
be verified in terms of eigenvalues and eigenvectors. In the framework 
of classical bifurcation theory, such an explicit characterization is very 
difficult, if not impossible, and the existence of heteroclinic orbits for 
a partial differential equation is often hard to prove. 

Although the focus in this article is on activator-inhibitor systems, 
the analysis is quite general and applies, for example, to systems con- 
sisting of a self-amplifying activator and a depleted substrate. 

Following is an outline of the paper. In Section [2], we formulate 
the reaction-diffusion problem for an activator-inhibitor mixture and 
rewrite it as an evolution equation in a function space. In Section [3l 
we study the exchange of stability, which is crucial for the stability 
and bifurcation analysis. The results of the bifurcation analysis are 
summarized for the one-dimensional case in Section H] and the two- 
dimensional case in Section [51 In Section [6l we illustrate the theoretical 
results on two examples, namely the Schnakenberg equation and the 
Gierer-Meinhardt equations. Section [7| summarizes our conclusions. 
Appendix |X] contains a brief summary of the attractor bifurcation the- 
ory from Ref. [11] and the reduction method introduced in Ref. [12j . 

2. Statement of the Problem 

Consider a mixture of two chemical species which simultaneously 
react and diffuse; one of the species is an activator, the other an in- 
hibitor of the chemical reaction. Their respective concentrations U and 
V satisfy a system of coupled nonhnear reaction-diffusion equations, 

Ut = fiU,V) + d,AU, 

^ ' ' Vt = giU,V) + d2AV, 
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subject to no-flux boundary conditions and given initial conditions. 
The functions / and g, which describe the kinetics of the chemical 
reaction, are generally nonlinear functions of the arguments. The dif- 
fusion coefficients di and ^2 are constant and positive. 

We assume that the system of Eqs. (12. ip admits a uniform steady- 
state solution which is positive throughout the domain. That is, there 
exist constants m > and v > Q such that 

(2.2) f{u,v)=^, 9{u,v) = ^. 

We are interested in solutions that bifurcate from this equilibrium so- 
lution and, in particular, in their long-term dynamics, under the as- 
sumption that the equilibrium solution (12.21) is stable in the absence of 
diffusion. 

For the bifurcation analysis, it is convenient to rescale time and space 
and rewrite the system (12.11) in the form 

Ut = ifiU, V) + AU, 
^ ' ^ Vt = ^g{U,V) + dAV, 

where 7 = 1/di and d = ^2/(^1. Thus, 7 is a measure of the ratio 
of the characteristic times for diffusion and chemical reaction, and d 
is the ratio of the diffusion coefficients of the two species. The above 
equations are satisfied on an open bounded domain, say Q C R" {n = 
1,2), while U and V satisfy Neumann (no-flux) boundary conditions 
on the boundary dQ of Q. 

2.1. Bifurcation Problem. Let 

(2.4) U = u + u, V = v + v. 

Since u and v satisfy the identities (12. 2p . we have 

f{U,V) = fu{u,v)u + f^{u,v)v + fi{u,v), 

g{U,V) = gu{u,v)u + g^,{u,v)v + gi{u,v), 

where /i and gi incorporate the higher-order terms in the Taylor ex- 
pansions. The functions u and v satisfy the equations 

ut = Au + -f{fu{u, v)u + f^{u, v)v) + 7/1 (m, v), 

vt = dAv + ^{gu{u, v)u + g^{u, v)v) + 75-1 (m, v). 

Henceforth we omit the arguments (m, v) and use the abbreviations 
for fu{u,v), et cetera. 

Since the variables U and V are associated with the activator and the 
inhibitor, respectively, of the chemical reaction, we have the inequalities 

(2.7) U > 0, g,< 0. 
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The equihbrium solution (12. 2p is stable in the absence of diffusion, so 
we also have the inequalities 

(2.8) fuQv - fvQu > 0, fu + Qv < 0. 

The first inequality in (12.81) . together with the inequalities (12. 7p . implies 

that f^Qu < 0. 

The problem as stated has two parameters, 7 and d. We represent 
the ordered pair by a single symbol, A = (7,(i), and consider A as 
the bifurcation parameter. The bifurcation is from the trivial solution, 
{u,v) = {0,0). 

2.2. Abstract Evolution Equation. The system of Eqs. (12. 6p de- 
fines an abstract evolution equation for a vector-valued function w : 

(2.9) ^ = L,w + G,{w),t>0; M^) = ( ; 

Here, Lx a linear operator in H of the form 

(2.10) Lx = -AD + -fB, 
where A : dom(A) ^ H is given by the expression 

-A 



f2.11) A = -AI 



-A 



on dom(A) = Hi = {w e {H\n)f : n ■ Vw = on dVl}, and B : H 
H and D : H H are represented by the constant matrices 

In Eq. (12. lip . A denotes the Laplacian, H'^{Q) is the usual Sobolev 
space, and the gradient on the boundary dfl of Q is taken component- 
wise. 

The nonlinear operator Gx : H —>■ H is given by 
(2.13) G,:^=(" ^7^^''"'"* 



Without loss of generality, we assume that Gx can be written as the 
sum of symmetric multilinear forms, 

(2.14) Gx{w) = Y,GxA^,--- ,^), 

k=2 

where Gk is a symmetric fc-linear form {k = 2,3,...). When the k 
arguments of Gk coincide, we write G^ with a single argument, Gkiw) = 
Gk{w, ... ,w). 
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The abstract evolution equation fl2.9p belongs to a class of equations 
analyzed in detail in Ref. [TT]; the relevant results are summarized for 
reference purposes in the Appendix, Section lA.li 



3. Exchange of Stability 
The inequalities fl2.8p imply that 

(3.1) det{B) > 0, tr(E) < 0. 

Under these conditions, diffusion has a destabilizing effect: At some 
critical value Aq of A, an exchange of stability occurs and the solution 
of Eq. fl2.9p bifurcates from the trivial solution. 

3.1. Eigenvalues and Eigenvectors of Lx and L*^. The negative 
Laplacian —A on a bounded domain Q G R" with Neumann boundary 
conditions is selfadjoint and positive in L^{Q). Its spectrum is discrete, 
consisting of eigenvalues pk with corresponding eigenvectors (fk, 

(3.2) -A(fk = pkipk, A; =1,2,..., 

We assume that the eigenvalues are ordered, < pi < p2 < ■ ■ ■ , and 
that the eigenvectors {fk}k form a basis in L^(f2). 

It follows from the definition (12.110 that A is selfadjoint and positive 
in H] its spectrum is also discrete, consisting of the same eigenvalues pk 
and the eigenvectors (pk once repeated. The operator Lx reduces via 
projection to its components on the linear span of each eigenvector of 
A. Let Ek be the component of Lx in the eigenspace associated with 
the eigenvalue pk, 

(3.3) EkiX) = -pkD + ^B, A; = 1,2,.... 

The determinant and trace of Ek{X) are 

, . det{Ek{X)) = 7Met(5) + 'ypklQvl - pkdinfu - pk), 
^ ■ ^ tr(Efe(A))=7tr(5)-pfe(l + d). 

Note that ii{Ek{\)) is negative everywhere in the first quadrant and 
becomes more negative as k increases. 

The eigenvalues of Ek{X) come in pairs, 
(3.5) 

PkiW = I (tr(Efc(A)) ± {{tT{Ek{X)f-Adet{Ek{X))f') , ^ = 1,2. 

They are either complex conjugate, with ^Pki = ^Pk2 < 0, or they are 
both real, with Pki + Pk2 < 0. We identify Pki with the upper (+) sign 
and Pk2 with the lower (— ) sign, so ^Pk2 < ^Pki- 
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(3.6) Wki 



The eigenvector corresponding to the eigenvalue Pki of Lx is 

ilfu - Pk- (3ki)^k 

The eigenvectors Wki and Wk2 are hnearly independent as long as I3ki 7^ 
(3k2- The set of eigenvectors {wki}k,i forms a basis for H. 

Note that B is not symmetric; its adjoint B* is the transpose B' 
of B. Hence, the adjoint of L\ is L\ = —AD + 'jB', and the adjoint 
of Ek{X) is -Efc(A) = —pkD + 7!?'. The eigenvalues of -^^(A) are f3ki, 
i = 1,2, the complex conjugates of the eigenvalues Pki of -E'fc(A) given 
in Eq. 03.51) . Since the latter are either complex conjugate or real, the 
eigenvalues of Lx and L\ coincide. The eigenvector corresponding to 
the eigenvalue Pki of L*^{X) is 



(3.7) w 



ki 



-19u'Pk_ 

ilfu - Pk- f3ki)'Pk 



3.2. Exchange of Stability. The equation det(ii^fc(A)) = defines a 

curve Afc in the (7, (i)-plane, 

(3.8) 

Ak = {(7, d) : det{Ek{\)) = 0} = {(7, rf) : d = 4(7)}, A; = 1, 2, ... , 
where 

(3.9) ,W) = 7M<'t(a)+-,p.|3„| 

Pkilfu - pk) 

The expression for dk can be recast in the form 

1n^ ^ / \ ^(s) Pk\fvgu\f (s)x_i , det(i?) {s)^-l 

(3.10) dk{l)-d^> = 75 — (7-7fc) + ^(7-7fc) , 

Ju Pkju 

where 7^''^ = pk/ fu and rf^**) = (det(5) + \fv9u\)/fu- This expression 
shows that (i) is symmetric with respect to the point (7^**'*, rf'^'*^) 
(ii) Afc has a vertical asymptote at 7 = pk/ fu] and (iii) A^ has an oblique 
asymptote with slope det{B)/{pkfu)- The symmetry point d'^'^^) 
is located in the first quadrant of the (7, d)-p\ane; 7^**^ increases as 
k increases, d^'^^ is independent of k. The vertical asymptote is in 
the right-half of the (7, (i)-plane, shifting to the right as k increases. 
The slope of the oblique asymptote is positive, decreasing to zero as 
A; increases. Therefore, each curve A^ has a branch in the positive 
quadrant of the (7, (i)-plane. The positive branches of the curves Ai 
and A2 are sketched in Fig. [H 

The curve A^ separates the region where det(i^^A;(A)) > (below the 
curve) from the region where det(£'fc(A)) < (above the curve). We 
focus on the region below the curve A2, bounded on the left by the first 
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Ri 

I 

Pi/fu 71 7 

Figure 1. Positive branches of Ai and A2. 



vertical asymptote at 7 = pi//„ and on the right by 7 = 71, where the 
curves Ai and A2 cross (see Fig. [1]). The curve Ai separates this region 
into two subregions, 

Ri = {\ = (7, d) : pi//„ < 7 < 71, < < c?i(7)}, 

i?2 = {A = (7, d) : pi/fu < 7 < 7i, c?i(7) <d< ^2(7)}- 

These regions are indicated in FigO 

Lemma 3.1. The eigenvalues (3ii (i = 1,2) of Lx satisfy the inequali- 
ties 



(3.12) 



3?/3ii(A) < 0, 3?/3i2(A) < ifXeRi, 
/3n(A)=0, /3i2(A)<0 z/AgAi, 
/5n(A)>0, A2(A)<0 ifXeR2. 



Furthermore, for k = 2,3, . . . , 
(3.13) 3fJ/3fci(A) < 0, 



k2 



(A) < ifXe Ai. 



Proof In we have tr(£'i(A)) < and det(Ei(A) > 0, so /?ii(A) and 
/?i2(A) are either complex conjugate with a negative real part, or they 
are both real and negative. On Ai, the leading eigenvalue /3ii(A) is 
zero. Since tr(i?i(A)) < 0, it must be the case that /3i2(A) is real and 



negative. In R2, we have tr(i?i(A)) < and det(£'i(A) < 0, so /3ii(A) 
and /3i2(A) are both real, and they have opposite signs. 

On Ai, det{E2{\)) > 0. Since det(£^fc(A)) increases with k, it follows 
that det(Efc(A)) > for A; = 2,3,.... Also, tr(Efc(A)) < 0. Hence, 
either /?fci(A) and /?fc2(A) are complex conjugate with a negative real 
part, or they are both real and negative. □ 

The lemma implies that all eigenmodes are stable as long as A is 
below the curve Ai. However, as soon as A crosses the "critical curve" 
Ai, the first unstable eigenmode appears and an exchange of stability 
occurs. 

4. Bifurcation Analysis - One-dimensional Domain 

We first consider the bifurcation problem (12.91) on a one-dimensional 
domain Q = (0, i). We reduce Eq. (12. 9p to its center- manifold represen- 
tation near a point Aq on the critical curve Ai, as proposed in Ref. [12] 
and sketched in the Appendix, Section IA.2I 

The eigenvalues and eigenvectors of the negative Laplacian subject 
to Neumann boundary conditions (see Eq. (13.21) ) are 



Pk = ^^(tt/^)^, v^fc(a;) = cos{k{7i/i)x), x G fi; A; = 1, 2, ... . 
The linear operator Lx decomposes into its components 



det{Ek{X)) = 7'det(5) + ^\g,\e{7i/if - dk^r^/ifi^U - k\n/if), 
tT{Ek{X)) = 7tr(5) - (1 + d)k\7:/ey. 



Each contributes two eigenvalues, (3ki and f3k2, to the spectrum 
of Lx] the expressions for P^i (i = 1,2) in terms of det(£'fc(A)) and 
tr(i?fc(A)) are given in Eq. (13.51) . The eigenvalues of the adjoint L\ are 
the complex conjugates, j3ki and J3k2- The eigenvectors of Lx and L\ 
corresponding to the eigenvalues [3ki and are 



(4.1) 



Ek{\) = -k\n/ifD + -fB, A; = 1,2,..., 



with 



(4.2) 



(4.3) 




(7/, 



u 



—7/1, COs{k{TT / i)x) 
k\7T/i)'^ - (]ki)cOs{k{7c/i)x) 



and 



(4.4) 




respectively. 
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4.1. Center-manifold Reduction. We are interested in solutions of 
Eq. (12.91) near a point Aq on the critical curve Ai. In the region Ri, just 
below Ai, both eigenvalues Pu and are real, with P12 < Pn < 0. As 
A approaches Aq, the leading eigenvalue Pn increases and, as A transits 
into i?2, Pii passes through and becomes positive. Thus, the first 
exchange of stability occurs. 

Lemma 4.1. Near the critical curve Ai, the solution of Eq. Ii2. 9\) can 
be expressed in the form 

00 

(4.5) w = yuWu + z, z = yuWu + XI XI y^i'^ki, 

k=2 i=l,2 

where the coefficient yn of the leading term satisfies the reduced bifur- 
cation equation, 

(4.6) ^ = p^^y^^ + ayl + oi\y^,\'). 

The coefficient a = a{X) are given explicitly in terms of the eigenfunc- 
tions of Lx and L\, 

(4.7) a(A) =a2(A)+a3(A), 
where 

< G'2(wii,^;2i),w*i >< ^2(^11),^;, > 

a3(A) = — < G3{wii),wl-^ > . 

Here, < ■ , ■ > denotes the inner product in H . (The subscript A on the 
k-linear forms has been omitted.) 

Proof. We look for a solution w of Eq. (12.91) of the form (14. 5p . In the 
space spanned by the eigenvector wn, Eq. (12. 9p reduces to 

(4.8) 

< > = < Lau;,w*i > + < Gx{w),wl^ > 



a2(A) = ^ V 



Pii < > yii + X < Gk{w),wl^ > . 



k=2 



To evaluate the contributions from the various terms in the sum, we 
use the asymptotic expression for the center-manifold function near Aq 
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given in the Appendix (Section IA.2p . Theorem I A. 2 1 
(4.9) 

v. = *i.(yu) = < ^''7''- ^ + o(l!,„l^), i = 2.3, . . . . 

(2/5ii - (3ki) < Wki, > 

The contribution from the bihnear form [k = 2) is 

< G2{w),wl^ > = < 6*2(1/11^11 + 2), > 
= < G2(wii),Wii > yf^ 

+ 2 < G2{wu,z),wl-^ > yu+ < G'2(z),<i > . 
The first term in the right member vanishes, because 

< G2(w^ii),Ki > = 0. 

The second and third term can be evaluated by means of the asymptotic 
expression (14.91) for the center manifold, 

< €2(1011, z),wl^ > = < €2(1011, W2i),wl^ > y2i + odyiiH 

i=l,2 

= |a2 < u^ii, w^n > 2/n + odyiiP), 
< G2(z),wl^ > = o(\yn\^), 

where 02 is defined in Eq. (14.71) . Putting it all together, we obtain the 
asymptotic result 

(4.10) < €2(10), wl^ > = a2< wu,wl^ > yl^ + o(\yii\^). 

The contribution from the trilinear form (A; = 3) is 

^^^^^ <G3(w),Wn >= <G3(wii),<i > yf^ + o(\yn\^) 

= a3<Wn,wl^>yf^ + o(\ynf), 

where as is defined in Eq. (14. 7p . The higher-order forms contribute 
only terms of o(|?/iip). □ 

4.2. Structure of the Bifurcated Attractor. The results of the bi- 
furcation analysis for one-dimensional spatial domains are summarized 
in the following theorem. 

Theorem 4.1. n = (0,i). 

• If a(Xo) < 0, then the following statements are true: 

(1) w = is a locally asymptotically stable equilibrium point of 
Eq. 113) for XeRi orXe Ai. 

(2) The solution of Eq. \2. 9\) bifurcates supercritically from 
(\q, 0) to an attractor A\ as A crosses Ai from Ri into R2. 
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(3) There exists an open set U\ d H with eU\ such that the 
bifurcated attractor A\ attracts U\\T in H , where T is the 
stable manifold of with codimension 1. 

(4) The attractor Ax consists of two steady-state points, 
and , 

(4.12) w7± = ±(/?n/|a|)i/Vi + ^A, A G i?2, 

where \\uJ\\\h = o{l3l(^). 

(5) There exists an e > and two disjoint open sets and 

in H, with G dU^ fl dU^ , such that G and 
linii^oo \ \w{t', Wq) — vj^Wh = for any solution w{t; wq) of 
Eq. I\2. 9\) satisfying the initial condition w{0]Wo) = Wq E 
and any X satisfying the condition dist(Ao, A) < e. 
• //a(Ao) > 0, then the solution of Eq. h2. 9\) bifurcates subcriti- 
cally from (Aq, 0) to exactly two repeller points as A crosses Ai 
from i?2 into Ri. 

Proof. Equation (14. 6 p shows that, if q;(Ao) < 0, then w = is a locally 
asymptotically stable equilibrium point. 

According to the attractor bifurcation theorem (Section IA.lt The- 
orem [XTT]), the system bifurcates at (Ao,0) to an attractor A\ as A 
transits from Ri into i?2- 

The structure of the attractor follows from the stationary form of 
Eq. (iD, 

All/11 + ai/ii + o(ll/iin = 0. 
The number and nature of the solutions of this equation does not 
change if the terms of o(|?/iip) are ignored, provided all solutions are 
regular at the origin. Thus, if a < 0, we find two solutions near |/ = 0, 

(4.13) y,, = ±{(3n/\a\f'^ + o{(3'l^). 

The last assertion of the theorem follows by time reversal. □ 

Theorem 14.11 shows that, if a(Ao) < 0, the attractor consists of two 
steady-state points, each with its own basin of attraction. The attractor 
bifurcation is shown schematically in Fig. [2j From the perspective of 
pattern formation, the theorem predicts the persistence of two types of 
patterns that differ only in phase; which of the two patterns is actually 
realized depends on the initial data. 

5. Bifurcation Analysis - Two-dimensional Domains 

Next, we consider the bifurcation problem (12.91) on a two-dimensional 
domain = (0,£i) x (O,^)- 
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Figure 2. One-dimensional domain: Supercritical bi- 
furcation to an attractor A = {w^,w^}. 

As in the one-dimensional case, we reduce Eq. (12. 9p to its center- 
manifold representation near a point Aq G Ai. 

The eigenvalues and eigenvectors of the negative Laplacian subject 
to Neumann boundary conditions (see Eqs. (14.31) and (14. 4p ) are 

y^kikiix) = COs{ki{TT/ii)xi)cOs{k2{lT/i2)x2), X = {Xi,X2) G Q. 

Here, ki and k2 range over all nonnegative integers such that \k\ = 
ki + k2 = 1,2,.... 

The eigenvalues I3k^k2i (^ = 1; 2) and the corresponding eigenvectors 
of L\ are given in Eqs. (13. 5p and (13.60 . respectively, where k now stands 
for the ordered pair (fci, k2). 

The dynamics depend on the relative size of £i and ^2■ On a rect- 
angular (non-square) domain, they are essentially the same as on a 
one-dimensional domain. For example, if £2 < ^1, then pio = (vr/^i)^ is 
the smallest eigenvalue of the negative Laplacian, with corresponding 
eigenvector ipio = cos((7r/£i)xi), and the leading eigenvalue of Lx is 
Pioi. This eigenvalue is simple, and the corresponding eigenvector is 

^c; n ,„ _f -7Acos((7r/£i)xi) A 

' ~ V ilfu - i^/^if - Aoi) cos((7r/4)xi) ) ■ 
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The center-manifold reduction leads to a one- dimensional dynamical 
system similar to Eq. fl4.6p . Lemma [4.11 and Theorem 14 . 1 1 apply verba- 
tim if Pii is replaced by /?ioi and Wu by Wioi everywhere. 

On the other hand, the dynamics become qualitatively different if 
the domain is square — that is, if = £2 = and fl = (0,£)^. The 
eigenvalues and eigenvectors of —A on the square are 

Pk^k, = {kl + kDin/if, 
<^kik2{x) = cos(A;i(7r/£)a;i) cos(/c2(7r/£)a;2), x = {xi,X2). 

Note that Pkik2 = Pfc2fci pair {ki,k2), so the eigenvalues Pkik2i 

{i = 1, 2) of L\ satisfy the same symmetry condition, Pkik2i = Pk2kii- 
To avoid notational complications, we consider two eigenvalues, even if 
they coincide because of symmetry, as distinct and associate with each 
its own eigenvector. Thus, we associate the eigenvector 

n\ _ / -lfv^k^k2 

^ ' " V - Pk.k2 - iik.k2^^k, 

with the eigenvalue Pkik2i, and the eigenvector 

^^'^"^ ~ ( {ifu - Pfcifca - hik2i)^k^k2 ) 

with the eigenvalue Pkik2i, whether ki and k2 are equal or not. 

5.1. Center-manifold Reduction. We are again interested in val- 
ues of A near the critical curve Ai, where the first exchange of stability 
occurs. The leading eigenvalues are /3ioi and Pqh. These eigenvalues 
coincide, but we consider them separately, each with its own eigenvec- 
tor. The two eigenvalues pass (together) through as A crosses Ai into 
i?2 from Ri, at the value A = Aq. 

Lemma 5.1. Near the critical curve Ai, the solution of Eq. Ii2. 9\) can 

be expressed in the form 

(5.4) w = yiWi + y2W2 + z, z = ^ ^ yk^k2iWk^k2h 

(fci,fc2):|A;|=2,3,...i=l,2 

where wi = wioi and W2 = woii- The coefficients yi and y2 of the 
leading terms satisfy a system of equations of the form 

Aoil/i + {ayl + cryDyi + o{\y\^), 

Pmm + {ayl + (Tyl)y2 + o(|y|^), 
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(5.5) 



dyi 
dt 

dm 

dt 



where /5ioo = Poii- The coefficients a = a(A) and a = cr(A) are given 
explicitly in terms of the eigenfunctions of Lx and L\, 

(5.6) a(A) =a2(A) +a3(A), 

where 



«2(A) = ^ V 



«3(A) = ^— — < Gs{wi),wl >, 

< Wi, wl > 

and 

(5.7) cr(A) = a2(A) + (T3(A), 

_ 4 ^ < G2{w2,Wiii),wl >< G'2(Wi,W2),Wii > 

<wi, wl > (2Aoi - Aii) < ^ui > 

3 

'^3(A) = < G3{wi,W2,W2),wl > . 

<wi,wl> 

The asymptotic estimates in Eq. / I5.5]) are valid as \y\ = ||yi || + ||y2|| 
0. 

Proof. We look for a solution w of Eq. (12.91) of the form (15. 4p . In the 
space spanned by the eigenvectors wi = wioi and W2 = wqu, Eq. (12.91) 
reduces to 

d °° 

< wi,wl > — = Pioi < wi,wl > yi + ^ < Gkiw),wl >, 

(5.8) 

d °° 

< W2,W2 > — = Pou < W2,W2 > y2 + ^ < Gk{w),wl > ■ 

k=2 



To evaluate the contributions from the various terms in the sums, we 
again use the asymptotic expression for the center-manifold function 
near Aq given in the Appendix (Section IA.2p . Theorem IA.21 

Ej=i,2 < G2{wj),wl^^^>y] 



(5.9) (2/3ioi - (3k,k2i) < Wk,k2i, wl,k2i > 



+ o{\y\'), h^k2, 



t \ < G2iwi,W2),wlk^> yiy2 , 2 

i^PlOl — Pkki) < Wkki,Wi^f,- > 

where = + ||l/2|P- 
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Consider the first of Eqs. fl5.8l) . Tlie contribution from tlie bilinear 
form is 

i=l,2 

+ 2J2 < G2iwi, z),wl> yi+ < G2{z), wl > . 

i=l,2 

The first term in the right member vanishes, because 

<G2iwi),wl> =0, i = l,2. 

The last term is asymptotically small, 

<G2{z),wl>=o{\yf). 

The second term involves an infinte sum over (fci, ^2) with \k\ = 2,3, . . .. 
Many of the coefficients are zero, because of the specific form of wi, 
W2, and Wk^k2j- The non-zero terms can be evaluated asymptotically 
by means of the expression fl5.9p . In fact, the only terms that are non- 
zero and contribute to the leading-order (cubic) terms in y are those 
with i = 1 and either {ki,k2) = (2,0) or (A;i,/c2) = (I,!)- Asymptotic 
expressions for 7/20* and ym (z = 1, 2) are given in Eq. fl5.9p . where we 
note that only the term with j = 1 contributes to y20i- 

Taken together, these observations show that the contribution from 
the bilinear form is 

(5.10) <G2{w),wl> = I <wi,wl> {a2yl + (r2yiy2)yi + o{\y\^), 

where 0^2 and a2 are defined in Eqs. (15. 6p and (15.70 . respectively. 
The contribution from the trilinear form is 

< Gs{w),wl >=J2< G3{wi),wl > + o(|?/n 

(5.11) i=l,2 

= < wi, wl > {a^yl + a^yDyi + o{\y\^), 

where as and are defined in Eq. (15. 6p and (15.70 . respectively. 

The computations for the second of Eqs. (15.80 are entirely similar. 
One finds the differential equation for y2 given in the statement of the 
lemma with the same expressions for the coefficients a and a. We omit 
the details. □ 

5.2. Structure of the Bifurcated Attractor. Before proceeding to 
the analysis of the structure of the bifurcated attractor, we recall the 
following result, the proof of which can be found in Ref. jl2] . 
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Lemma 5.2. Let yx G &e a solution of the evolution equation 

^ = Xy-G,Ay) + o{\y\''), 

where G\^k is a symmetric k-linear field, k odd and k > 3, satisfying 
the inequalities 

Ci\y\'^'<<Gx,kiy),y><C2\y\'^' 

for some constants C2 > Ci > 0, uniformly in A. Then yx bifurcates 
from {y,X) = (0,0) to an attractor Ax which is homeomorphic to . 
Morover, one and only one of the following statements is true: 

(1) Ax is a periodic orbit; 

(2) Ax consists of an infinite number of singular points; 

(3) Ax contains at most 2{k + 1) singular points, which are either 
saddle points or (possibly degenerate) stable nodes or singular 
points with index zero. The number of saddle points is equal to 
the number of stable nodes, and both are even (2N, say). If the 
number of singular points is more than 4N (4N + n say, where 
AN + n < 2{k + 1)), then the number of singular points with 
index zero is n and N + n > 1. 

The results of the bifurcation analysis for two-dimensional spatial 
domains are summarized in the following theorem. 

Theorem 5.1. VL = {0,i)^. 

• If a{\o) < and a(Ao) + cr{\o) < 0, the following statements 
are true: 

(1) w = is a locally asymptotically stable equilibrium point of 
Eq. ^ for XeRiorXe Ai. 

(2) The solution of Eq. \2.9) bifurcates from (Aq, 0) to an at- 
tractor A{X) as X crosses Ai from Ri into R2. 

(3) The attractor A{X) is homeomorphic to . 

• If a{Xo) < and cr(Ao) < 0, the attractor A{X) consists of an 
infinite number of steady-state points. 

• If a(Ao) < and cr(Ao) > 0, the attractor A{X) consists of 
exactly eight steady-state points, which can be expressed as 

(5.12) wx = Wx + uJx, AGi?2, 

where Wx belongs to the eigenspace corresponding to /?ioi and 

W^xWh = o{\\Wx\\h)- 

• If a{Xo) > and a(Ao) + cr(Ao) > 0, the solutions of Eq. h2. 9\) 
bifurcate from (Ao,0) to a repeller A{X) as X transits into R2. 
Also, A{X) is homeomorphic to S^. 
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Proof. Equation flS.Sp shows that, if a(Ao) < and a(Ao) + cr(Ao) < 0, 
then 1(7 = is a locally asymptotically stable equilibrium point. 

It follows from Lemma [5l2] and the attractor bifurcation theorem lA.lt 
that the system bifurcates at (Aq, 0) to an attractor A\ as A transits 
from Ri into R2, and that Ax is homeomorphic to S^. 

The structure of the bifurcated attractor is found from the stationary 
form of Eq. (15.51) . Ignoring the terms of o(|?/p), we have the system of 
equations 

(Aoi + ayl + (^vDvi = 0, 

(5.13) 2 2 

(/3ioi + + (Jy^)y2 = 0. 

If a < 0, a + 0" < 0, and a > 0, the system (I5.13P admits eight nonzero 
solutions near y = 0, 

2/1=0, yj = f3^o^/\a\; 

(5.14) 2/2 = 0, yl = p^oi/\a\; 

yl = yl = Aoi/I" + 

These solutions are regular, so Eq. (15. 5p also has eight steady-state 
solutions; they differ from the solutions of Eq. (15.131) by terms that are 

The last part of the theorem follows by reversing time. □ 

Theorem 15.11 shows that, if a(Ao) < and a(Ao) + o"(Ao) < 0, the 
bifurcation is an S^-attractor bifurcation. If both a(Ao) and cr(Ao) are 
negative, the attractor consists of an infinite number of steady-state 
points; on the other hand, if a(Ao) < and cr(Ao) > 0, the attractor 
consists of precisely eight steady-state points. Figure [3] shows the phase 
diagram on the center manifold after bifurcation, when A has crossed 
the critical curve Ai into the region R2. The phase diagram consists of 
eight steady-state points and the heteroclinic orbits connecting them. 
The odd-indexed points (Pi, P3, P5, and P7) are minimal attractors; 
they correspond to striped patterns. The even- indexed points (P2, P4, 
Pe and Pg) are saddle points. 



6. Examples 

We illustrate the preceding results with two examples from the the- 
ory of pattern formation in complex biological structures, namely the 
Schnakenberg equation ^19j and the Gierer-Meinhardt equation [5] ; see 
also Ref. P [TOl [131 [151 [ISl [13 [23] . 
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Figure 3. Two-dimensional domain: ^^-bifurcation 
with eight regular steady states. 

6.1. Schnakenberg Equation. A classic model in biological pattern 
formation is due to Schnakenberg [19] . 

Ut = j{a-U + U'V) + AU, 
Vt = -fib- U^V) + d/\V, 



on an open bounded set C R*^ (n = 1, 2), with Neumann boundary 
conditions and given initial conditions. The constants a and h are 
positive; 7 and d are positive parameters. The system admits a uniform 
steady state, 



a + h 

b 



(6.2) 

(a+6)2 

The Schnakenberg equation is of the type fl2.9p . with 

and a nonlinear term Gx of the form fl2.13p . with 

^) = 7 — TTx?"^ + 2(a + h)uv + u'^v, 
[a + by 

gi{u,v) = -fi{u,v). 
The conditions (12. 7p and (12. Sp are satisfied if 

a<6, 6-a<(a-h bf. 
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6.1.1. One-dimensional Domain. ^ = (0, 1). 

An evaluation of the inner products in Eq. (14. 7p with the MAPLE 
software package yields the expression 

M^ S1 + S2 + S3 
aiX) = ^, A e Ai, 

-2726(a + 6) + (7^ - pi) 

where 

Si = h'(7 + Pi){a + &)'((5pi + P2i)ia + b) + 7(2a - b)) 

^ (7(2a -b)+ 2pi(a + &))(4pi + 7 + ^ ^ ^ 

/32.(2726(a + 6)-(7^-p2-/?2.n ' ' ' ' 

S3 = §7^(7 + Pi)(a + bf (7(& - a) - pi(a + b)) . 

Note that a(A) depends only on 7 if A G Ai; we use the short-hand 
notation 0(7) = q;(7, ^1(7)). 

6.1.2. Two-dimensional Domain. Q = (0, 1)^. 

An evaluation of the inner products in Eqs. (15. 6p and (15. 7p with the 
MAPLE software package yields the expressions 

aiX) = 7T, A G Ai, 

-2726(a + &) + (7fef-poi)' 



and 



a(X) = r , A G Ai, 

^' -2726(a + 6) + (7^-poi)2' 



where 

s^ = |7'(7 + Pi)(a + &)'((5pi + f32o^)ia + b) + 7(2a - b)) 

/32o.(27^&(« + &)-(7^-4pio-/520iP)' ' ' ' 
S3 = 77^(7 + Pi){a + b)^{-f{b -a)-pi{a + b)); 



s 



+ 5)4(^ + p^)((2a - 6)7 + {(3,,, + 3pi)(a + b)) 
^ ((2a - &)7 + 2(a + &)pi)(7 + 2pi + Pm) ^ 

X — ■ ■ — 7 rrrr 1 ^ — J-i-^i 



= §7^(7 + Poi)(a + &)^(7(& -a)- Poi(a + &))• 

6.1.3. Numerical Results. Numerical results are given for a = |, 6 = | 
in Fig. m and for a = 2 and b = 100 in Fig. [51 In the former case, 
there is no bifurcation; in the latter, there is a pitchfork bifurcation at 
Ao = (7o,c?(7o)). 
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Figure 4. Schnakenberg equation with a — |, b — |. 

(i) Positive branches of Ai and A2 in one and two di- 
mensions; (ii) Graph of a in the one-dimensional case; 
(iii) Graph of a in the two-dimensional case; (iv) Graph 
of (7 in the two-dimensional case. 




Figure 5. Schnakenberg equation with a = 2, b = 100. 
(i) Positive branches of Ai and A2 in one and two di- 
mensions; (ii) Graph of a in the one-dimensional case; 
(iii) Graph of a in the two-dimensional case; (iv) Graph 
of (7 in the two-dimensional case. 
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6.2. Gierer— Meinhardt Equation. Another model for the forma- 
tion of Turing patterns was proposed by Gierer and Meinhardt [5j, 



(6-3) ,.,2 



Ut = 7(a -bU + UyV) + AU, 
Vt = 7(f/2 ~V)+ dAV, 



on an open bounded set Q C R" (n = 1, 2), with Neumann boundary 
conditions and given initial conditions. The constants a and b are 
positive. The system admits a steady state, 



(6.4) 



g+l 



b2 

The Gierer-Meinhardt system is an equation of the type (12. 9p with 

/ (l-a)& b2 

D — I 1+a (l+a)2 
I 2(l+a) _^ 

The nonlinear term Gx is obtained by expanding around the steady- 
state solution, 

/i M, = u - ——uv 

(1 + ay \ 1 + a 

(1 + a)^ \ 1 + a 

giiu.v) = u^. 
The conditions (12. 7p and (12.81) are satisfied if 

a < 1, 6 < . 

1 — a 



6.2.1. One- dimensional Domain. f2 = (0, 1). 

An evaluation of the inner products in Eq. ( 14. 7p with the MAPLE 
software package yields the expression 

S1 + S2 + S3 
a\A) = 2! 

-1'^. + (7+? - p: 
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where 

^\2a - 1)62 + 76(p, + 2a{h^ + 5pi)) + 2pi(l + a){P2^ + 4pi) 



= - ^7^^ 



S3 = 2 



(l + a)8 

^\2a - 1)9 + 7K4pi(l + a) + /?2^) + 2p^ (l + a) 
3 7^65(^(1 _ a)fe - pi(l + a))(7a6 + pi(l + a))^ 



i = 1,2 



5 '^•) 



(1 + a) 



2 ' 



6.2.2. Two-dimensional Domain. Vt = (0, 1)^. 

An evaluation of the inner products in Eqs. fl5.6p and (15. 7p with the 
MAPLE software package yields the expressions 

a{X) = 

and 

s^ + s'^ + s^ 
a{X) = 25 

-1'^, + (7^ - fH«] 



where 

s. = - h'b' 



2 

-f%^ + 7&(pio + 2a{/320i + 5pio)) + 2pio(l + a){P20i + 4pio) 

(l + a)8 

72(2a - 1)62 ^ ^b{Ap^o{l + a) + /?2o^) + 2p?o(l + a) 



S3 



/?20. (-7^i^+(7^-Pio-/32o.)') 
3 7255(7(1 - a)6 - pio(l + a)){^ab + pio(l + a)f _ 



^ = 1,2, 



(1 + a) 



^2(2a - 1)62 ^ ^^(^^^ + 2a(/3in + 3pio)) + 2pio(l + a)(Ai, + 2p 



10 J 



(l + a)8 

^2(2a _ 1)52 + ^^(^^^^ + 2pio(2a + 1)) + 2p2o(l + a) 



Pu. (-7^T^+(7^-pio-/3n.)') 
3 _ o7'&^(7(l - - Pio(l + a)){^ab + pio(l + a)f 



1 = 1,2, 



=3- 



(l + a)8 
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6.2.3. Numerical Results. Numerical results are given for a = |, 6 = 1 
(Fig. [6]) and a = | and & = | (Fig. [7]). In the former case, there is no 
bifurcation; in the latter, there is a bifurcation at Aq = (70,^(70)). 




Figure 6. Gierer-Meinhardt equation with a = ^, b = 
1. (i) Positive branches of Ai and A2 in one and two 
dimensions; (ii) Graph of a in the one-dimensional case; 
(iii) Graph of a in the two-dimensional case; (iv) Graph 
of a in the two-dimensional case. 



7. Conclusions 

In this paper we considered the evolution of an activator-inhibitor 
system consisting of two morphogens on a bounded domain subject to 
no-flux boundary conditions. Assuming that the system admits a uni- 
form steady-state solution, which is stable in the absence of diffusion, 
we focused on solutions that bifurcate from this uniform steady-state 
solution. The bifurcation parameter A represented both 7, the ratio 
of the characteristic times for chemical reaction and diffusion, and d, 
the ratio of the diffusion coefficients of the two competing species (ac- 
tivator and inhibitor). We showed that, for such a system, there exists 
a critical curve Ai in parameter space such that, as A crosses Ai, a 
bifurcation occurs (Lemma 13.11) . 

While a linear analysis around the uniform steady state suffices to 
obtain information about the formation of patterns, a nonlinear anal- 
ysis is needed to gain insight into the long-time asymptotic behavior of 
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Figure 7. Gierer-Meinhardt equation with a = |, 6 = 
|. (i) Positive branches of Ai and A2 in one and two 
dimensions; (ii) Graph of a in the one-dimensional case; 
(iii) Graph of a in the two-dimensional case; (iv) Graph 
of a in the two-dimensional case. 

the solutions after bifurcation. This issue is intimately connected with 
the long-term persistence of patterns. 

In this paper we used the theory of attractor bifurcation, in com- 
bination with a center-manifold reduction, to analyze the long-time 
dynamics of bifurcated solutions. We considered two cases: diffusion 
on a (bounded) interval or a (non-square) rectangle, and diffusion on a 
square domain. In the former case, we showed that a bifurcation occurs 
as A crosses a critical curve Ai, and the bifurcation is a pitchfork bi- 
furcation. Theorem 14.11 gives an explicit condition for the existence of 
an attractor. The attractor consists of exactly two steady-state points, 
each with its own basin of attraction. The two steady states correspond 
to patterns that differ only in phase; which of them is eventually real- 
ized depends on the initial conditions. Essentially the same conclusion 
holds in the case of diffusion on a rectangular (that is, non-square) do- 
main; in particular, roll patterns emerge as a result of the bifurcation. 

In the case of diffusion on a square domain, the dynamics are qualita- 
tively different. Theorem 15.11 gives explicit conditions for the existence 
of an S'^-bifurcation. The bifurcated object consists of either an infinite 
number of steady states, or exactly eight regular steady-state points 
with heteroclinic orbits connecting them. Thus, two types of patterns 
may arise; for example, in the formation of animal coat patterns, we 
might expect stripe patterns or spot patterns. 

25 



Thus, in both the one- and two-dimensional case we have given a 
complete characterization of the bifurcated attractor and, therefore, of 
the long-time asymptotic dynamics of the bifurcated objects. 

Appendix A. Attractor Bifurcation and Reduction 

Methods 

In this appendix, we summarize the attractor bifurcation theory of 
Ref. [TT] and the reduction methods introduced in Ref. [T2|. The func- 
tional framework is that of two Hilbert spaces, H and Hi, where Hi is 
dense in H and the inclusion Hi "-^ H is compact. 

A.l. Attractor Bifurcation Theorem. Let A : Hi ^ H he a linear 
homeomorphism, and let L\ : Hi ^ H he a, compact perturbation of A 
which depends continuously on a real parameter A, 

(A.l) La = -A + 5a. 

The operator Lx is sectorial; it generates an analytic semigroup S\{t) = 
{e^^^}t>o- Fractional powers are defined for all a G [0, 1]; the do- 
main of is dom(L^) = Ha, where Hq = H and Ha^ C Ha^ if 
a2 < ai- 

Let Gx : Ha ^ H he a nonlinear C'-bounded map (r > 1) for 
some a G [0,1), which depends continuously on A and satisfies the 
asymptotic estimate 

(A.2) Gx{w) = o{\\w\\hJ, 

as \w\ — 0, uniformly in A. 

Let Wx{t; Wq) G if be a solution of the initial value problem 

(A.3) ^ = LxW + Gxiw),t>0; w(0) = u^o, 

at 

where wo & H is given. In terms of Sx, we have 

(A.4) wx{t]Wo) = Sx{t)wo+ [ Sx{t-s)Gx{wx{s;wo))ds, t > 0. 

Jo 

Definition A.l. (i) A set G H is a (positive) invariant set of 
Eq. (Hl^j if Sx{t)Ti = E for any t > 0. (ii) An invariant set G H of 
Eq. ^A.!^} is an attractor z/S is compact and there exists a neighborhood 
U G H ofTj such that 

lim diist Hiwxit; Wq), E) = 

t— >oo 

for any wq G U. (Hi) The largest open set U satisfying the above 
condition is the basin of attraction o/S. 
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Definition A. 2. (i) A solution {wx,X) of Eq. liA.3\) bifurcates from 
(0, Ao) if there exists a sequence of invariant sets {fin} of Eq. 
with ^ Qn such that 

lim max \ w\ = 0, lim dist(An, Aq) = 0. 

n— >oo ui£f2n n— >oo 



(a) If the invariant sets Qn o'^e attractors of Eq. then the bifur- 



cation is called an attractor bifurcation. (Hi) If the invariant sets Qn 
are attractors and are homotopy equivalent to an m- dimensional sphere 
, then the bifurcation is called an S^-attractor bifurcation. 

Tlie theory of Ref. [TT] is developed under the conditions that (i) the 
spectrum of A is discrete and consists of positive eigenvalues pk with 
corresponding eigenvectors ipk & Hi, 

(A.5) ALpk = Pk^k, fc = l,2, 

where < pi < p2 < ■ ■ ■ ^ oo; (ii) the vectors {(pk}k form an orthog- 
onal basis of H; and (iii) there exists a constant 6 G (0, 1) such that 
Bx '■ Hg —>■ H is bounded uniformly in A. 

The eigenvalues of Lx are /3fc(A), k = 1,2,.... The assumption is 
that there is a critical curve that separates the region Rm where 
(the real parts of) the first m eigenvalues are negative from the region 
Rm+i where (the real parts of) the first m eigenvalues are positive, and 
an exchange of stanility occurs if A = Aq G A^. 

Theorem A.l (Attractor Bifurcation Theorem fllj). Let l3k{X), k = 
1, 2, ... , be the eigenvalues of Lx (counting multiplicity). Suppose that, 
at some critical value A = Aq, the first m eigenvalues /3i(A) through 
/9m(A) cross the imaginary axis into the right half of the complex plane, 
while the remaining eigenvalues, l3k{X) for k = m + l, m + 2, . . remain 
in the left half of the complex plane. Let Ei be the eigenspace of Lx at 

m 

Ei=[j{weHi: {Lx, - l3k{\o)yw = 0, ^ = 1, 2, . . .} , 

fe=i 

and let w = be a locally asymptotically stable equilibrium point of 
Eq. / CO) at \ = Xq. Then 

(1) Eq. liA.3\) bifurcates from {w, A) = (0, Aq) to an attractor Ax as 
A transits from Ri into R2, dim^A ^ [m — l,m], and Ax is 
connected if m > 1; 

(2) Ax is the limit of a sequence of nested m- dimensional annuli 
Mi with Mj_|_i C Mi; in particular, if Ax is a finite simplicial 
complex, then Ax has the homotopy type of S"^"^ ; 
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(3) For any wx € Ax, Wx can be expressed as 

wx = Wx + zx, WxeE^, zx = o{\\Wx\\Hy, 

(4) There is an open set U C if with E U such that Ax attracts 
U\T, where T is the stable manifold ofw = with codimension 
m. 

A. 2. Reduction Method. A useful tool in the study of bifurcation 
problems is the reduction of the equation to its local center manifold. 
The idea is to project the equation to a finite-dimensional space after 
a change of basis; details can be found in [12j. 

Let A be close to a critical value Aq. Suppose that the spaces H and 
Hi are decomposed, 

(A.6) H = Ei® E2, Hi = Ei® E2, 

where Ei and E2 are invariant subspaces of Lx, Ei is finite dimensional, 
El = El, and E2 is the closure of E2 in H. The decomposition flA.6P 
reduces Lx, 

(A.7) Lx = Cx,i © Cx,2, 

where £a,i = Lx\ei : Ei Ei and Cx,2 = Lx\e2 ■ E2 E2. If the 
decomposition is such that the real parts of the eigenvalues of Cx,i 
are nonnegative at A = Aq, while those of Cx,i are negative, then the 
solution Wx of Eq. flA.3l) can be written as 

(A.8) w = W + z, WeEi,zeE2, 

where W and z satisfy the system of equations 

dW 

-— = Cx,iW + gx,iiW,z), 
(A.9) ^1 

= C-x,2Z + gx,2{W,z), 

with Qx,i = PiGx, Pi : H ^ Ei being the canonical projection. 

By the classical center-manifold theorem (see, for example, Refs. O 
[23]). there exist, for all A sufficiently close to Aq, a neighborhood Ux C 
£^1 of = and a center-manifold function : Ux ^ Ei, which 
depends continuously on A, such that the dynamics of Eq. flA.3p are 
described completely by the dynamics of the finite-dimensional system 

^ = Cx,iW + GxAW, ^\W)), WeUxCEi. 
dt 

The following theorem gives an asymptotic approximation for as 
A — >^ Aq or, alternatively, as — >■ 0. The proof of the theorem is given 
in Ref. [I2]. 
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Theorem A. 2. Assume that G\, the nonlinear part of Eq. liA.3\) . 
IS C°°, Gx{w) = YlT=p'^\k{w) for some p > 2, where Gx,k{w) = 
G\^k{w, . . . ,w), and G\^k : Hi x . . . x Hi ^ H is a k-linear map. Then, 
under the conditions of Theorem \A.l\ the center-manifold function $^ 
can he expressed as 

+ o{\\W\n, We El, A^Ao, 
where j3 = {pi, ... ,13^). 
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